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Abstract 

The nonet symmetry scheme seems to describe rather weU the masses and 77—7/ mixing angle of the ground 
state pseudo-scalar mesons. It is expected that nonet symmetry should also be valid for the matrix elements 
of the pseudo-scalar density operators which play an important role in charmless two-body B decays with rj 
or rj' in the final state. Starting from the divergences of the SU{2>) octet and singlet axial vector currents, we 
show that nonet symmetry for the pseudo-scalar mass term implies nonet symmetry for the pseudo-scalar 
density operators. In this nonet symmetry scheme, we find that the branching ratio B PP, PV, with 77 in 
the final state agrees well with data, while those with rj' are underestimated, but by increasing the B rj' 
form factor by 40 — 50%, one could explain the tree-dominated B~ tt't]' and B^ — > p^rj' measured 
branching ratios. With this increased form factor and with only a moderate annihilation contribution, we 
are able to obtain 62 x 10^^ for the penguin-dominated B^ K^-q' branching ratios, quite close to the 
measured value. This supports the predicted value for the B rj' form factor in PQCD and light-cone sum 
rules approach. A possible increase by 15% of (0|s?75s|ss) for 7/0 would bring the predicted B^ K^rj' 
branching ratio to 69.375 x 10^^, very close to experiment. 

PACS numbers: 13.25Hw 
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I. INTRODUCTION 



Unlike the low- lying vector mesons where the flavor diagonal 1~ qq states are eigenstate because 
of the OZI selection rule, the 0~ pseudo-scalar qq state can mix with each other. Since QCD 
interactions through the exchange of gluons are flavor- independent, one expects the wave function 
for the pseudo-scalar nonet also flavor-independent in the limit of vanishsing current quark mass 
{nig 0,q = u,d, s) and the r] and rj' can be described as two linear combinations of the qq state, 
the SU(3) singlet r/o and the SU(3) octet 773 which mix with each other through a small 5*^7(3) 
symmetry breaking mixing parameter. In fact, with m„ and <C rUg , rUg <^ Aqcd , and because 
of the U{1) QCD-anomaly, the 770 mass is much larger compared to the r/g mass, the rj — rj' mixing 
angle is 0{ms/ Aqcd) so that the physical r] and rj' are almost pure ryg and r/o eigenstate respectively, 
in contrast with the ideal mixing for the 1~ low-lying vector meson states. Another features of 
the 0~ qq nonet is that, because of the spontaneous breakdown of SU{3) x SU{3) symmetry, the 
octet mesons are massless Goldstone bosons in the limit of vanishing current quark mass. This 
simplifies considerably the description of the ground state pseudo-scalar meson system. As shown 
in[l|, a rather accurate description of the mass and mixing angle in the rj — rj' system is obtained 
by adding an U{\) QCD-anomaly term for the 770 rnass in the nonet pseudo-scalar mass matrix. 
This mass matrix is generated by the quark mass term and is the leading term in the large 
expansion while higher order terms in the chiral Lagrangian|^ is 0{\/N(.) and is thus suppressed in 
the large N^. limit. This justifies the nonet symmetry mass term for the pseudo-scalar mass matrix. 
Vice- versa, from the nonet symmetry value for the off-diagonal mass term < rjQ\H^^\ri^ >, where 
^.e = ™. + ™„ r» + ™, U one wo.ld .et . mixi.. an* . = - 18" i„ good a^eemen. with a 
value e ^ -(22±3)° in [l], or 6* « -(18.4±2)° in j3| and a similar value 9 « -(17-20)° [4] obtained 
from the two-photon decay width of rj and rj' . However, if we use the Gell-Mann-Okubo(GMO) 
mass formula for the octet mass rrig, we would have, in terms of the -q — rj' mixing angle 6 

rn^ = m| — tan O'^ {irif^i — rrig) (1) 

which, for 6 = —18° gives = 483 MeV, about 60MeV below experiment. Thus the rj — rj' 
mixing which contributes to L7 in f2i has driven the m„ below the GMO value by 63 MeV. This 

nn 

is also the case with a nonet mass matrix in the quark basis [5|, l6[ which has a large rj — rj' mixing 
and an upper bound for the rj mass far below experiment. The higher order terms L4, L^, Lq, 
and chiral logarithms obtained in Ref. [3] shift upward by a similar amount with the result 
that the rj mass is very close to the GMO value, in agreement with experiment. Similar result is 
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also obtained in [6] more recently. Thus, nonet symmetry seems to be a good approximation for 
the 0~ nonet mass term. One could then go further and try to see if the matrix elements of the 
pseudo-scalar density local operator e.g. si^^ s could also satisfies nonet symmetry. This will allow 
a simple calculation of the penguin matrix elements in the charmless two-body decays of B meson 
with r] or rj' in the final states. In this paper we will use the divergence equation for the octet and 
singlet axial vector current to show that nonet symmetry scheme for the mass term implies nonet 
symmetry for the pseudo-scalar density qi^^q for ij and 7]'. The basic idea is to include in the 
matrix elements of the axial vector current and its divergence the r/o,8 pole contribution which will 
add the mixing mass term < 7?o|-f^SB|??8 > to the divergence equation and allows us to obtain the 
nonet symmetry expression for the matrix element of the pseudo-scalar density operators between 
the vacuum and 779,8 • In the next section we will first derive a divergence equation for the u 7^75 u 
and 57^^755 axial vector current, in the presence of the SU{3) x S'C/(3)-breaking Hsb current 
quark mass term. Section III is an analysis of Pr] and B^ Pi]' ,P = K^,Tr^ in QCD 

Factorization (QCDF) with nonet symmetry for the pseudo-scalar density and B ^ rj and B ^ rj' 
transition form factors. We find that the branching ratio for modes with rj in the final state agrees 
well with data, while those with rj' in the final state are underestimated. We then increase the 
B ^ rj' form factor by 40 — 50%, to bring the tree-dominated B~ ir^r]' and B^ p^ij' to the 
measured values. The increased form factor is then used to obtain a branching ratio close to data 
for the penguin-dominated B K~r]' decay. 

II. PSEUDO-SCALAR DENSITY MATRIX ELEMENT AND NONET SYMMETRY 

Let \r]o >, |r/8 > be the SU(3) singlet and octet eigenstate of the / = 0, pseudo-scalar nonet in 
the absence of the SU{3) symmetry breaking quark mass term Hsb, in terms of the flavor diagonal 
qq component: 

\r]Q >= {\uu + dd + ss >)/a/3, 

\r]s>= {\uu + dd-2ss>)/V6. (2) 

Consider now the matrix element of the axial vector current matrix element m7^75U and S'j^'j^s 
between the vacuum and r/o and r/8 : 

< 0|u7^75u|r/o >= ifuPf,/V3, 

< 0|n7^75n|r/8 >= ^fuP^J./VQ■ (3) 
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and 

< 0|s7^75s|?7o >= ifsPfi/Vs, 

< 0\s-f^j5s\r]s >= -2i /sP^/Vq. (4) 

where fu and fs are defined as the decay constant of uu and ss state respectively. Except for the 
momentum dependence factor p^, the above axial vector current matrix elements depend on the 
same L and fs according to nonet symmetry scheme with identical qq spatial wave function in ryo 

M 

and %[1], but fs could be different from fu by an SU{ci) breaking s-quark mass term. The octet 
Ag^ and singlet j4o^ axial vector current matrix elements between the vacuum and 7?8i% are then 
given by: 

^n\A \ ^ ifu + fd + ^fs) 

<0\A^o\m>= ^ Pa*- (5) 

(p^ is the 4-momentum of rjQ and r/g . Similar, for other members of the SU{3) octet, we have 
/tt and fx for 7r+ = ud, = us meson respectively. Assuming each s-quark contributes to the 
decay constant a symmetry breaking term e, to first order in e, (Rewriting fqq = /g), we haveQ]: 

fK = fus = (1 + e) fud , 

/, = (1 + 2e) + e)/K. (6) 

The usual way to obtain the pseudo-scalar density matrix elements is to take the divergence of the 
axial vector current between the vacuum and the pseudo-scalar meson octet. For example, taking 
the matrix elements of uij^d, uij^s, [uij^u — (dij^d) between the vacuum and vr"*", K^, i:^, 
respectively, we have: 

fnBo{mu + md) = {■mu + md){0\uij5d\ud), 

fKBo{mu + rus) = {ruu + ms){0\ui'y5s\us) . (7) 

and for tt^ : 

fuBo{mu + rrid) = {niu + md){0\uij5u\uu) . (8) 



where t 



le TT and K meson masses are the usual expressions in terms of Bq and the current quark 



mass 



[ll, 2I]. The expression for ir^ is obtained by putting: [s 



(0|iZi75'u|7r") = -(0|(ii75d|7r'^). (9) 



Apart from the difference in f-^ and fx, we see that the above pseudo-scalar density matrix element 
in Eq.([7]) and Eq.dS]) satisfies SU(3) symmetry. We will see below that to have nonet symmetry 
for the pseudo-scalar density matrix element between the vacuum and 7/0,85 the pole term in the 
divergence equation must be included. We now consider the divergence of the I = A^^j, and 
axial vector current: 

^sm = s7m75S- (10) 

The divergence is given by: 

— a ~ 

dA^ = 2(m„m75U -|- niddi'y^d) + 2—GG, 

4tt 

OAs = 2mssij5S + G. (11) 

The matrix elements of dA^ and dAg between the vacuum and r/0,8 ^-re given by: 

(0|5An|%>=2m„(0|ni75n|r/o) + 2md{0\d ij^dlrjo) , (12) 
+ 2{0\^GG\m)- 

(y. ~ 

(0|5^|r?o)=2m,(0|si75s|??o) + {0\-^G G\r]o). (13) 

47r 

and for r/g : 

(0|5^n|??8>=2m„(0|ni75n|r/8) + 2mrf(0|(Ji75(i|r78>, 

+ 2(0|^GG'|r?8) , (14) 

(0|9^|r?8)=2m,(0|si75s|??8) + {0\^GG\r]s). (15) 

47r 

In the limit niu = rrid = , since the l.h.s of Eq. (|14|) is /um|, the matrix element 2{0\^G G\7]s) 
on the r.h.s is 0(mg) and is given by the r/o pole contribution. We now evaluate Eq. (jl2til3|) and 
Eq. (jl4m5p with the pole terms included using the nonet symmetry expressions for rriQ g and rngg [ll 
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mj = Bo- {2ms +m), 

ml = ttiq -h -Bo^(ms 2m), 
2 

"^08 = BQ-V2{-ms + m). (16) 

in standard notation]^ (m = (m„ + md)/2). At the r/o and rj^ mass, = ttiq and = m| in the 
l.h.s of Eq. (|12til3|) and Eq. ()14m5|) respectively. As mentioned above, since m^^ <C mg-, SU{3) is 



broken and the 7?o,8 pole will contribute to both the l.h.s and r.h.s of Eq. (jl2H13|) and Eg. (j 1411 15]) 

. The pole terms on the r.h.s come from the QCD-anomaly matrix element {0\^GG\7]q) and 

{0\j^G G\r]s) induced by 5C/(3)-breaking r/o — rnixing mass term rngg . The presence of the 

r/0,8 pole term is important, since its contribution is the same order as the current-quark mass 

terms in ijiq g . Indeed had we dropped the r/o,8 pole term we would run into contradiction with 

the divergence equation. To obtain the pseudo-scalar density matrix elements, let us bring the 

p^-dependence pole term in the l.h.s to the r.h.s of Eq. (jl2til3]) and Eq. p^ - (fT5|) . Putting = fd 

and (0|ui75'U,|r?o,8) = {0\di'y5d\r]0,8) , we find, for tjq : 

12 11 2\/2 1 

fu^ifhl + Bo-{ms + 2m)) = /„-^?fio -/«-^(-Bo-^(m - m^)) 2--^m(0|u i75'u|uu),(17) 

12 12 2\^^ 1 

^''^^^o + ^o-^i^s + 2rh)) = -/^— 5o-^(m - m^) + 2— m5(0|si75s|ss). (18) 

and. similarly, for rjs : 

12 1 2\/2 1 

fu-j={Bo-{2ms + m)) = -fu-^Bo—^{m-ms)+2—m{0\ui-f5u\uu), (19) 

2 2 1 2V2 2 _ _ 

-fs^{Bo-{2ms + m)) = -fs-j=Bo^-{m - tUs) - 2-j=ms{0\si-f5s\ss). (20) 

Comparing the l.h.s and the r.h.s of Eq. (fT7|) and Eq. ifTH]) . we get the pseudo-scalar density 
matrix element for ryo: 

{0\ui-f5u\uu) = Bofu, (21) 
{0\sij5s\ss) = BoU (22) 

Similarly, by comparing l.h.s and the r.h.s of Eq. ()19p and Eq. (l20p . we get the same expression for 
the pseudo-scalar density matrix element, but in ryg . 

We have shown that, by including the ijq and r]s pole in the divergence equations, and by 
using the nonet symmetry expressions for the current quark mass contributions to the ijq and 
ijs mass, the pseudo-scalar density operators matrix elements between ijq and r/g can be ob- 
tained by nonet symmetry and quark counting rule. Like the matrix elements {0\uij5d\TT^) , 
(0|n i75n|7r'^) and (0172 275517^"^), they are given by the parameter Bq and the decay constant in- 
volved. Experimentally, from the known value of the ij — r]' mixing angle, 6 = (—20 it 2)° , one has 
TTigg = —(0.81 lb 0.05) to be compared with the nonet symmetry value of uiq^ ~ — 0.90m^[l|, 
we expect nonet symmetry for the pseudo-scalar density matrix elements in rj — ij' valid to this 
accuracy. Since the octet m| mass gets about 15% increase from higher order terms L4, L5, Lq, Lg 
and chiral logarithms^], Eqs. (|19ti20]) show that (0|si75s|ss) in 77 will be increased by a similar 



amount. Note that the r.h.s of Eqs. (jl9M20l) gets this increase from higher order terms in the pole 
and other terms. Higher order SU{3) breaking contribution to the singlet mass is not known, 
but if we assume a similar 15% increase from the nonet value in Eq. (jl6p . {0\si^5s\ss) in r/o will 
also be increased by a similar amount. This could be another source of enhancement for the 
B Krj' branching ratio, as found below. We note that it might be possible to obtain the pseudo- 
scalar density matrix elements in Eqs. (|17tiT8]) and Eqs. (|19ll20|) using the known values of ttiq g and 
rriQg , but because of the precise dependence on quark mass is not known and experimental errors 
involved, the physical interpretation of the result will be lost. We would like to stress that in our 
derivation, the anomaly contribution to the ??o mass has been included in the divergence equa- 
tion, thus the enhancement factor for (0|s i75s|r/o) suggested in |9i] would have the origin elsewhere. 
With the pseudo-scalar density matrix elements given above and nonet symmetry for the B ^ ri,r]' 
transition form factors, we shall now compute the B^ — > K^rj, K^r]' and B^ — > tt^t/, vr^r/' decay 
branching ratios in QCD Factorization(QCDF). 



III. B- if (77,77 ) AND B- Tr-{f],i) DECAY IN QCD FACTORIZATION 



The B M1M2 decay amplitude in QCDF is given by[lO|, [ill]: 
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A{B ^ M1M2) = -^Yl ^pbK ^ - E 4{MiM2\Oi\B)H + fBfMjM.bi , (23) 

^ p=u,c \ i=l i J 

where the QCD coefficients contain the vertex corrections, penguin corrections, and hard specta- 
tor scattering contributions, the hadronic matrix elements {MiM2\0i\B) h of the tree and penguin 
operators Oi are given by factorization model[l^. [l^ , h are annihilation contributions. The values 
for a^,p = u,c , computed from the expressions in [10) [ill 



at the renormalization scale 



nib, 



with rrih = 4.2 GeV are: 



''10 



*10 



-0.033 - 0. 013 i + 0.0009 pH, 

-0.027 -0.017i + 0.0009pH, 

-0.045 - 0.003 i, = -0.042 - 0.013 i, 

-0.0004 - 0.0001 i, = 0.0004 - 0.0001 i, 

-0.0011 -O.OOOH -0.0006 

-0.0011 + 0.0006 i - 0.0006 



(24) 
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for i = 4,6, 8, 10. For other coefficients, af = = ai : 

ai = 1.02 + 0.015* - O.OUpH, 

02 = 0.156 - 0.089i + 0.074^^, 

03 = 0.0025 + 0.0030 i - 0.0024 pH, 
05 = -0.0016 - 0.0034i + 0.0029 pn, 
aj = -0.00003 - 0.00004 i - 0.00003 
09 = -0.009 - 0.0001 i + 0.0001 pn- 

(25) 

where the complex parameter pn exp(i0//) represents the end-point singularity contribution in the 
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ll| (we have put the phase (pn = 



hard-scattering corrections Xh = (1 + pn ewii'pH)) ^^i^j^) 
in the above expressions). 

For the annihilation terms, we have: 

62 = -0.0038 - 0.0065 pA - 0.0018^0^, 

63 = -0.0065 - 0.0150 PA - 0.0085 p% 

bf" = -0.00011 - 0.00015 PA +0.000003 p^. (26) 

where bi are evaluated with the factor fsfMifivh included and pA , like pH, appears in the divergent 
annihilation term Xa = (1 + pAexp(i0yi)) ln(^^). 

For the CKM matrix elements, since the inclusive and exclusive data on \Vub\ differ by a large 
amount and the higher inclusive data exceeds the unitarity limit for Rf, = \Vudyub\/\^<^^cb\ '^ith 



data. We havejl5t|: 



the current value sm(2,) ^ 0.687 . 0.032^, we shah determine W the more prec.se ,^.1 



Vub I = I sin /3 W 1 + . (27) 

\V^,\ V sm^a 



ud 

With a = (99t^^)°Q and \ Vcb\ = (41.78 ± 0.30 ± 0.08) x lO'^ Q, we find 



iKfel = 3.60 X 10"^ (28) 

in good agreement with the exclusive data in the range \Vub\ = 3.33 - 3.51Q . The recent 
measurements of the Bs — Bg mixing also allow the extraction of \Vtd\ from Bd — Bd mixing 
data. The current determination [17] gives |Vtd/Vts| = (0-208^Q'QQg) which in turn can be used to 
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determined the angle 7 from the unitarity relation [l^: 



\^cbV*a\ I. , cos^a 



\V*, 



Vtd\ = sin7A/l + (29) 



;le a = 91.8°, in good agreement with the value 



with \ Vti,\ = 1, we find 7 = 66° which implies an ane 
found in the current UT-fit value of (88 it 16)° [181]. In the following in our B decay calculations, 
we shall use the unitarity triangle values for \Vub\ and 7. For other hadronic parameters we use 



the values in Table 1 of [ill and take m^(2GeV) = 80MeV, /„ = f^, f, = f^(l + 2(^ - 1) 



■/vr 

For the i? — > vr and B ^ K transition form factor, we use the current light-cone sum rules central 
value Q : 

F(f^(0) = 0.258, Fo^^(O) = 0.33 (30) 
With r] — 7]' mixing angle 9 = —20°, we have 

\r]) = (0.58(|nn) + \dd)) - 0.57|ss)), 

\r]') = {OAO{\uu) + \dd)) + 0.82|ss)). (31) 



F(f^(0) = 0.58Fo^^(0), F^'^'iO) = 0.40Fo^^(0). (32) 



From Eq.(|3H we find : 



The B K{r]',r}) decay amplitude can now be obtained from the factorization formula for the 
hadronic matrix elements in Eq. (j23p with the pseudo-scalar density matrix element obtained in 
Eq. (j22p and the form factors given above. We have 

{0\s i-f5s\rj) = CriBofs, 

{0\sij5s\r]') = C\>Bofs. (33) 

where Bq = m\/{ms + m) and C^^ = —0.57, C^' = 0.82, the fraction of ss state in rj and rj' 
respectively. This contributes to the Og matrix element a term fs^^, with = 2m|^/(mb + 
rh) {rus + rh), similar to that in B^ — > K^ir^ decay, except that in B~ — > K~r] and B^ K~r]' , 
the Og matrix element is enhanced by a factor fs/fx- In this way, the decay amplitude in unit of 
GeV are : 

AiB'^TT^TT^) = (0.110 + 0.204?) X lO"'^, 
A{B^^K-TT+) = -(0.368 + 0.004i) X 10~^(Fo^^^(0)/0.258) 

-(0.090 + 0.002 i) X lO"'^. (34) 



from which the branching ratios are, with pn = 0,PA = 0.6 (only the central values for the relevant 
parameters are used in the calculations) 

B{B- vr^vr") = 5.050 x 10~^ 

B{B^ K-TT+) = 18.249 x 10^^ (35) 
in good agreement with the current measured branching ratios [3] 

13{B- 7r-7r°) = (5.7 ± 0.4 x 10"^ 

B{B° K-TT+) = (19.04 ± 0.6) x 10"^ (36) 



We note a sizable annihilation contribution, given by the last term in Eq. ()34p . is needed to produce 
a large B{B^ K^ir^) . This is not surprising since annihilation contribution is also needed to 



explain the large branching ratios of B^ tt^K*^ and B^ ^ K decay 2l|]. Our result also 
shows that the values 0.258 for F^'^{0) and 0.33 for F(f^(0) given above are reasonable. We will 
use these values in the calculation of the decay modes with r] and r]' . We find 

A{B- -^K-ri) = -(0.283 + 0.032 i) x IQ-"^ (F(f ^''(0)/0.150) 
+(0.317 + 0.080i) X 10"'^(Fo^^^(0)/0.33) 
+(0.015 + 0.0004 i) X 10"^. (37) 

A{B- K'T]') = -(0.192 + 0.022 i) x 10"^ (F(f (0)/0.104) 
-(0.425 + 0.039 i) x 10"^ (Fo^"*^(0)/0.33) 
-(0.111 + 0.003 i) X lO"'^. (38) 

where the last term in Eq. (j37p and Eq. (j38p are the annihilation contributions (/)// = 0,pA = 0.6). 
The predicted branching ratios are then : 

B{B- K-7]) = 0.431 X 10-^ 

BiB-^ K-7]') = 48.263 x 10"^ (39) 
to be compared with the current experimental values [2^: 

B{B- K~rj) = (2.2 ± 0.3) x 10-^ 

B{B- K-T]') = (69.7l|?) X 10"^ (40) 

We see that the B{B~ — > K~r]') is underestimated by about 30%, while the B{B~ K~r]) is 
very much suppressed, but because of large cancellation in the B~ — > K~r] amplitude due to the 
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negative ss amplitude in the r] meson wave function, a precise prediction for B{B~ K~rj) is 
more difficult. For K~r]', since 63 contributes both to — > K~~Tr^ and B^ K^rj' 

decays, it is difficult to adjust the annihilation term for B^ — > K~r]' without overestimating the 
B^ K^ir^ branching ratio. Another possibility is to increase the form factor Fq~*^ (0) from 
the nonet symmetry value to bring the predicted value closer to data. That this is the case can be 
seen by looking at the B~ — > vr^r/' decays. We have: 

A{B-^TT'r]) = (0.119 + 0.147 i) X 10"' (F(f^''(0)/0.150) 
-(0.002 - 0.003 i) X 10"^ (Fo^^''(0)/0.258) 
-(0.004 - 0.003 X lO"'^. (41) 



A{B- -^TT-ri') = (0.081 + 0.100 x lO"'^ (F(f ^'''(0)/0.104) 
+(0.008 - 0.002 i) X lO"'^ (Fo^^''(0)/0.258) 
+(0.033 - 0.021 i) X lO"'^. (42) 

(the last term in the above amplitudes is the annihilation contributions). This gives: 

B{B- ir-T]) = 3.388 x 10"^ 

B{B- TT-7]') = 1.910 X 10"^ (43) 
comparing with the current measured branching ratios [2^: 

B{B~ TT~r]) = (4.4 ± 0.4) x 10~^ 

BiB- ^ TT-ri') = (2.6l°|) X 10"^ (44) 

we see that the predicted B{B^ t^ v) agrees more or less with experiment, considering theoretical 
uncertainties in the CKM parameters and in the B ^ tt and B ^ K form factors, while B{B 



■K'T]') is below the Babar value of (4.0 ± 0.8 ± 0.4) x 10"*^ [20(. Existing QCDF calculations [23] 
also underestimate B{B~ — > p~r]') by a factor of ~ 2 as seen from the recent data|2j] which gives: 

B{B- p-7]) = (5.4 ± 1.2) X 10-^ 

B(B- ^ p-r/') = (g.ltli) X 10"'- (45) 

Since the above tree-dominated decays with rj, rj' in the final state are more sensitive to the F^^^ 
and F^^^ form factor, by increasing the F^^"^ form factor by 40 — 50% from the nonet symmetry 
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value, one could bring ^ vr rj'), B{B p r]'), and B{B K r/'), closer to the measured 

branching ratios. For example, by taking F^^^ (0) = 0.156, one gets: 

B{B- TT'T]') = 3.888 X 10-^ 

BiB~ K'T]') = 61.837 x 10^^ (46) 

which largely improves the prediction for B{B~ — > K~rj') but the predicted B{B~ — > iT~ri') slightly 
exceeds the HFAG new average, though consistent with the Babar value for this mode. We note also 
the predicted B{B^ P v') 2^ approaches the measured value with the increased form factor 
Fq~*'^ (0) = 0.156 . As mentioned earlier, additional source of enhancement of B{B^ K^t]') 
could come from a possible higher order SU{3) breaking effects in the matrix element {0\s i^5s\ss) 
for r]Q. Assuming a 15% increase of this matrix element from its nonet value, we would have 
B{B^ K^T]') = 69.375 x 10~^, very close to the measured value. 

IV. CONCLUSION 

We have shown that nonet symmetry for the pseudo-scalar meson mass term implies nonet 
symmetry for the pseudo-scalar density matrix element. We then use nonet symmetry for the 
pseudo-scalar density matrix element and the i? — > ry, i? — > 77' form factors to compute two- 
body charmless B decays with rj, rj' in the final state. The discrepancy with experiment for tree- 
dominated decays with rj' in the final state indicates that the F^^^ (0) form factors should be bigger 
than the nonet symmetry value by 40 — 50%. This value together with a moderate annihilation 
contribution found in B^ K^ir^ decay, produces B~ K'rj' branching ratio close to data. 
Our value for the Fn^^ (0) form factor supports the current calculations in PQCD and light-cone 



sum rules approach 



24l |. A possible increase by 15% of the pseudo-scalar density matrix element 



(0|s i75,s|,ss) for r]Q would bring the predicted B K ri' branching ratio very close to experiment. 
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